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ABSTRACT 


Mao lasoeeoumatalytical and numerical solutions to the 
meopwem Ol inviscid-jet deflection from plane and axisym- 
metric concave surfaces are presented. The methods of solu- 
erom fOr the two-dimensional cases included a novel method 
developed in the present investigation, as well as those 
due to Levi~Civita and Riemann and Hilbert. These methods 
have provided indirect solutions for the shape of the curved 
solid boundary in terms of a given jet~departure and jet- 
deflection angle. The axisymmetric cases, namely, the deflec- 
Tion of axisymmetric jets from hemispherical thrust reversers 
are solved directly through the use of the finite element 
method and a novel iteration scheme. The real-fluid effects 
on the jet deflection are estimated by comparing the results 
with those obtained experimentally. 

The analytical and numerical methods developed or used 
herein are sufficiently general to yield direct and indirect 
Solucvitons for one general nozzle and curved-reverser combi- 
nations of Special interest to thrust reversal on aircraft 


engines. 
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PN RODUCTION 


Micmomce me oroiw ot tie detleclion of finite jets by 
Eermetral Or curved deflectors is a matter of importance 
mimtcme GCesign of target—type thrust reversers on jet-aircraft 
engines, of flip-buckets on spillways, of Pelton wheels, etc. 
In general, the deflection of a free jet by a solid boundary 
is well suited to potential-flow analysis because of the 
dominance of inertia and pressure intensity in the establish- 
ment of the flow pattern. The design of impulse machinery, 
thrust reversers, etc. utilizing this momentum change could 
be facilitated greatly if the idealized geometry of the system 
under potential flow conditions were known because such condi- 
tions represent asymptotic values which are approached as the 
effects of secondary variables such as entrainment, boundary 
layer, compressibility, jet attachment to adjacent surfaces, 
ever ale Cdeereased. With such information available, refine= 
ments of design could be based upon a secure knowledge of the 
fundamentals, and many rules of thumb could be replaced with 
Preeecise Guanricarive data. Specifically, if the total angle 
through which the jet is deflected is determined for conditions 
GiebOun partial and complete interception by the boundary, then 
the principle of impulse and momentum can be used to compute 
forces or other dynamic characteristics of the system. 

The’ problems discussed herein fall in the general category 


of finite or infinite, two-dimensional or three-dimensional, 


Ik 





Tse ouomer Aanviscid,; and compressible or incompressible flow 
mesure DOd@es= of arbitrary shape. In fact, in its broadest 
sense, they constitute the most fundamental and also most 
eomplex problems of fluid dynamics. Cur inability to deter- 
mine Ueemcerparauron pelmnvus, or lines, on a body and to estab-— 
imooea link between #he separation and the unsteady fluid 
motion in the wake of the body gave rise to numerous exact 
and approximate methods of analysis. Most of these analyses 
rely in part or whole on the experimentally observed charac- 
teristics of flow to determine the magnitude of the disposable 
Perametvers introduced into the analysis. Suffice it to say 
oiao future efforts will continue to be based upon partly a 
secure knowledge of the fundamentals of fluid mechanics and 
partly on experimental facts. 

me rdemely, the introduction of simplifying assumptions such 
esmemvescid, incompressible flulds, the use of bodies with 
a salient edge at which the separation point is fixed and 
minaity the restriction of the analysis to either two-dimen- 
sional or axisymmetric flows enable one to obtain solutions 
which are often of practical significance and in conformity 
with the experimental observations. 

The present work deals with inviscid, incompressible flows 
of finite extent impinging upon two-dimensional or axisymmetric 
bodies with sharp edges, and makes use of all the available 
exact or approximate methods of analysis in arriving at novel 
solutions towards the prediction of the characteristics of 


deflected jets. 
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A closer look at both the problem at hand and the methods 
available requires the further specification of the geometri- 
Ga Comm ifurations Of both the nozzle from which the jet is 
emerging and the body upon which the jet is impinging. The 
classical theory of jet deflection deals essentially with 
nozzles and bodies composed of straight segments. This is 
primarily because the succesSive transformations used in the 
analysis lead to relatively integrable equations for the 
segmental boundaries and to exceedingly complex integral 
equations in the case of curved boundaries. Furthermore, 
the fact that there cannot be a corresponding complex-function 
theory for three-dimensional space, the limitations of our 
Gmgersctanding of the characteristics of axisymmetric or, in 
general, three-dimensional jet-deflection problems becomes 
guite clear. 

imorder tO bring into closer focus the results of the 
past efforts as well as the limitations of the existing 
Metniogds, we Will, in the foregoing, describe briefly the 
Suggestions made concerning the "rounding of the corners" 
of solid boundaries, describe the essential features of the 
problems encountered through the use of a specific example, 
and then discuss the history of the two-dimensional jet 
deflection problem. The comments concerning the previous 
work on axisymmetric jet deflection will be reserved for 
Section III where such cases are studied through the use of 


the finite element method. 





Peer BRIEF REVIEW OF SEGMENTAL JET DEFLECTORS IN 
TWO~DIMENSIONAL FLOWS 


The two-dimensional counterpart of the jet-deflection 
meoolem Nas been treated by several investigators through the 
use of the powerful analytic-function theory and successive 
conformal transformations. Sarpkaya [1] solved the U-shaped, 
two-segpment, deflector problem where the turning angle between 
the segments is limited to -90 degrees. Tinney et al. [2] 
extended this analysis to the case where the turning angle 
between the symmetrically situated segments is preater than 
90 degrees. Later, Chang and Conley [3] presented an analysis 
for a bucket composed of a series of segments of arbitrary 
number, length, and angles; however, the basic as well as 
Beoecvsacal problem of the direct analysis of jet deflection 
by curved buckets remains unsolved primarily because of the 
Meaenhematvical difficulties encountered. It is partly for the 
piempese Of L1llustrating these ditficulties and partly for 
the purpose of familiarizing the reader with the fundamental 
concepts and the analytical transformations that the basic 
steps in the analysis of the deflection of a two-dimensional 
jet from a U-shaped segmental deflector will be discussed. 

The flow in the physical z-plane (see Fig. 1) may be 
transformed into a hodograph-plane through the use of 

-i6 dw ae: 


BN ele = aS 


The flow in the hodograph plane may be transformed by means 
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(a) z-plane 


(b) ¢-plane 





(c) Q2-plane 
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(d) t-plane 


Figure 1 Transformation planes for the U-shaped bucket 
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of Planck's logarithmic transformation into an Q-plane 


where { is defined by, 
+ id (2) 


The Q-plane is composed of straight lines. The horizontal 
jines represent the solid boundaries and the vertical lines 
the free surfaces. Such a polygon may be transformed (e.g. 
see Milne-Tomson [4]) through the use of the Schwarz-Chris- 
toffel transformation onto either the upper or the lower 
half of a t-plane by writing 


QO = M, {$+ ap No (3) 


ai a5 Tt 
(t-t,) (t-t.) 


iimee won the t values in the t-plene may be chosen arbitrarily 
in accordance with the Schwarz-Christoffel transformation and 
the rest must be assigned unknown parameters. Accordingly, 
the points E, A, C are assigned the values -k, k and f 


respectively. Then, Eq. (3) becomes 


Ct Gat (4) 


Re Nf eee Oh 
Gye 2 


After integration and evaluation of the constants M,and N, 


one obtains 


b@ie i 2) and. ta. =k 


(>) 


——_— 





and 


Q = dn —— + for -k< t<l (6) 


the potential function W in the upper half of the t-—plane 

may be written in terms of the sources and sinks representing 
wae ilow. Noting that there is a source at D and a sink at 

C and that only the upper half of the t-plane is used, one 
has 


V.d 
W = — [ln (t-f) - ln (t+t1)] (7) 


meom this relationship and the definition of t in Eq. (1) 
and 2 in Eq. (2) one obtains the final expression that relates 


Che z-plane to the t-plane, 


psec a Pe 


Zz -{ Vi-t? tN we [ ae 1} AG (8) 
d [a> Co ee 


‘N 


Since the purpose of this section was only to show the general 
peoccaure msed for this type of UCransformavtions, the results 
of the integration are not shown here. The details are given 
by Sarpkaya [1]. 

It is apparent from the foregoing that whenever the solid 
boundaries in the physical plane (see Fig. 1) are composed of 
straight segments, then the N-plane also is composed of 


straight lines. In other words, the flow retains a constant 
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direction along the solid boundaries (lines parallel to the 
logarithmic axis in the Q-plane), and a constant velocity 
along the free surfaces (lines normal to the logarithmic 
axis in the Q-plane). Since such a polygonal boundary can 
always be transformed onto a t-plane through the use of the 
Schwarz-Christoffel transformation, all two-dimensional jet- 
deflection probiems of this nature may be, at least theoreti- 
meby, solved. 

It is also apparent from the foregoing that whenever both 
the magnitude and the direction of the velocity vary along a 
meena boundary, Chen the 2—-plane is not in general composed 
of straight lines or of lines conveniently intersecting each 
other at 90 or 180 degree angle. In general if the curved 
Memmaaries in the physical plane is given, then the corres- 
ponding portion of the Q-plane will be an unknown curve 
since the magnitude of the velocity is not initially known. 
Miemmay, therefore, attempt to simplify the problem, for 
example, by specifying the variation of the magnitude of the 
velocity along the boundary. This does not in any way alle- 
viate Mice rculbres tO be encountered later since vnere 
does not exist a transformation function comparable to that 
paver DY ochwarz and Christoffel which could transform an 
Q-plane composed in part by curved segments onto a t-plane 
in which the potential-function may be written. Consequently, 
one will either seek other methods of handling the curved 
boundaries or completely abandon the direct approach of 


Seraaning a solution for a given geometry. The indirect 
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m@aemoach Called into action will seek a family of solutions 
for a family of initially unknown curved boundaries and then 
will let the designer choose one among those obtained to suit 
meeoeGesired needs. In the following, the efforts made by 
mers if rounding the corners of buckets otherwise composed 
Memoeroiehit Jines and the efforts to obtain indirect solutions 


mo be described. 


B. ROUNDING OF THE CORNERS AND BOUNDARIES 

Because of its theoretical as well as technological im- 
meiemace, the problem of rounding the sharp corners on rigid 
boundaries over or through which the fluid flows has attracted 
The attention of many mathematicians and hydrodynamicists. 
mesor tO a brief exposition of these efforts, it should be 
noved that there has not been, up to now, a sufficiently 
general method of studying the flow around curved bodies and 
Wat the techniques so far developed are not only difficult 
to handle but also require the use of special techniques for 
eeeeieven problem. 

1. Curved Walls of Constant Pressure 

As cited earlier, both the magnitude and the direc- 

tion of velocity vary along a curved boundary and the func- 
tional relationship between the wall geometry and the vector 
Pewocivy 1s not known a priori. However, by making a simpli- 
fying assumption such that the pressure or the total velocity 
along the curved wall remains constant, one can round the 


corners in a bucket otherwise composed of straight segments. 


any, 





The method may best be described by applying it to the U- 
wmeoead bucket case Giscussed earlier. 

Let the Q2~plane (see Fig. lc) be modified in such a 
manner so that it takes the shape shown in Fig. 2a. Evidently 
the magnitude of velocity along PQ remains constant and the 
Girection varies from t/2 at P to m7 at Q. The Q-plane is 
still composed of straight lines and may, therefore, be 
Bralstormed a2nto either the upper or the lower half of a 
t-plane through the use of the Schwarz-Christoffel transfor- 
mation. The remainder of the analysis will not be presented 
here. Suffice it to say that the transformation now involves 
pyo additional terms in the form of Vt-t,, Vt-t, and, 
therefore, renders the succeeding integrations relatively 
more difficult to perform. Be that as it may, the method 
is in principle capable of generating curved corners along 
which the pressure or the velocity remains constant. There 
Zomoutirtelent flexibility in the analysis to move the points 
A and @ along the boundaries to obtain smaller or larger 
muevecmocorners Within reasonable limits. 

| mecurved Walls with Special Pressure Distributions 

A bolder approach to the modification of the Q-plane 
is its modification in such a manner that while the direction 
of flow along the curved wall varies linearly, the magnitude 
of the velocity varies logarithmically. This procedure still 
allows the “%-plane to be composed of straight segments (see 


Fig. 2b). It is immediately apparent that the application 
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eps: -plane for 
constant pressure -in(q/Vj) 
curvature 





D 
(b) Q-plane for 
variable pressure =I (q/ Vay) 
auEVvabure 





feiee- or Z2-plane (d) t-plane with modifications at P and ¢ 


Figure 2 Modifications of the 2- and t-planes to round the corne: 
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of the Schwarz-Christoffel transformation now yields terms 
such as 


0. TI O/T 
fiat.) (t-t,) 


as a consequence of the fact that the line segments in the 
f—-plane, no longer intersect each other normally. Only 
Mmcer rare circumstances, i.e. for special values of as 
Ong Sat that the integral equation resulting from the Schwarz- 
Christoffel transformation may be integrated in terms of 
the Beta-functions. In general, the transformation may be 
antegrated only numerically. Since the determination of 
the coordinates of the rigid boundaries in the physical 
plane requires another integration, similar to Eq. (8), the 
solution of the problem via a two-stage numerical integration 
Mecomes not only difficult but also quite approximate. 

This method has been applied by Schieldrop [5] to 
the analysis of Borda-mouth-piece type flows with curved 
edges. The method is not in general easy to apply and requires 
fmeeueay Ceal Of ingenuity and sophistication in numerical 
integrations. 

Seeeeurved Walls Through the Modification of the t—plane 

Cockcroft [6] introduced a method by means of which 
the corners of a solid boundary may be slightly rounded by 
modifying the t-plane. Consider the flow in a 90 degree 
corner as shown in Fig. 2e. This plane may be regarded as 


either z or the N-plane. In the absence of round corners, 


fees 





Phe Schwarz—-Christoffel transformation relates this plane to 


the t-plane (see Fig. 2d) by 


Zorg = 1 [Meo ae + N, (9) 
(t-a) Vi¢=b 


The corner at B may be rounded in an indirect way by replacing 
{to in the above equation by a a AY t-1 aha tas 
Peon rcuLion the argument of Vt which changed tori? 2uat 

t = Q now varies smoothly when t moves from t = -1 to t = tl. 
ma@meouner words, a curve is produced in the Zz or N=-plane whose 


curvature depends on the arbitrary value of ir. It should be 


noted that for a value of A, say 


ene velocity or the pressure along the curved wall remains 
Setiscanmr. In fact, evaluating dW/dz in the range of —-l<t<1, 


one obtains 








This method was successfully used by Cockcroft [6] 
in calculating the electrical-stress distribution around 
certain two-dimensional conductors. A brief discussion of 
mas method and the difficulties associated with its appli- 


cations are also described by Carrier et al. [7]. 
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y ewevine OL the Boundaries Through the Use of 
pocetal Transformations 


The Riemann—-Hilbert transformation has been exten- 
mm ely used by Larock and Street [8] in the analysis of 
flow about supercavitating hydrofoils and by Larock [9] in 
eaalyzZing in an indirect manner the efflux from curved 
m@eezles. the mathematical details of this method will be 
Mmememesca further later in connection with its application 


memcie jeu Geflection problem. 


Peeealyols Of Flow about Curved Bodies by 
Levi-Civita's Method 


Because of its challenging nature, the question of 
the existence and uniqueness of potential flows past general 
obstacles, having boundaries of given shape, has intrigued 
many outstanding mathematicians such as Brillouin [10], 
Villat [11], Levi-Civita [12], and others. In fact, Levi- 
Civita may be said to have solved the inverse problem of 
fescribing the class of all jets divided by curved barriers. 
meetouel No specific solutions have been presented, the prob- 
tem of determining such flows was reduced in special cases 
to the solution of nonlinear integral equations with appro- 
Pelate boundary conditions. The solutions are obtained in 
general by solving such equations numerically and by suitably 
choosing the necessary number of parameters. Such parameters 
described for example the variation of the inclination of 
the velocity vector, the variation of the magnitude of veloc- 
ieyesetc. it is important to know that the present problem 


of the deflection of a two-dimensional jet from a curved 
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Boundary has not been previously solved by any of the methods 
memerar described. 

ite oculonwlitLevi=Civata's method, among others, 
fmmebe Gaescribed in greater detail in generating a family 
Seecurved buckets for a given jet and prescribed jet 


@erebkection angle. 


C. SUMMARY 

Major emphasis in the foregoing introduction has been 
Omeune direct or indirect analysis of the impingement of 
ideal fluids on two-dimensional straight or curved boundaries. 
These analyses are often difficult to apply and the amount 
of work required may not be commensurate with the need for 
eesolution. In such cases, the use of relatively less- 
m@mecty Mevnods may be adequate or even necessary. For two- 
dimensional flows in general and for axisymmetric flows in 
particular approximate methods have been developed. A de- 
Talled discussion of these methods with the exception of the 
finite element method may be found in Robertson [13]. These 
methods may be classified as graphical (field plot or flow 
net), finite-difference numerical analysis (iteration, relax- 
ation, Monte Carlo), and analogies (electrical analogy, Hele 
Shaw, and membrane analogy), and finally the finite element 
method. The first three of these approximate methods have 
been widely used and amply described by Robertson [13]. The 
finite element method which came into existence during the 


past decade is currently being applied to the solution of 
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two-dimensional and axisymmetric Laplace-field problems. 
Peoetbailed discussion of this method and its application to 
axisymmetric jet deflection problems will be presented in 
meevton Jil. In the following, exact methods of analysis 


of the deflection of two-dimensional jets from curved buckets 


will be undertaken. 


26 





IIT. ANALYSIS OF TWO-DIMENSIONAL JET 
DEFLECTION FROM CURVED BUCKETS 


iveune following, three methods will be employed to 
obtain a family of curved—bucket shapes in terms of a given 
Pomeccettecction angie and the angle of departure at the lip of 
the two-dimensional bucket. The jet will be assumed to be 
free and impinging symmetrically upon the bucket. These 
methods consist of the modified hodograph-method, Levi- 
mevowa Ss method, and the Riemann-Hilbert method. The first 
of these is novel and the other two have been previously 


employed as discussed in Section I.B. 


A. MODIFIED HODOGRAPH METHOD 

The method consists of the definition of the hodograph- 
mere. 1ts transformation into an infinite strip, and subse- 
quently, into a t-plane through the use of the Schwarz-Chris- 
rorirel transformation. The flow in the t-plane is expressed 
iaeverms of appropriate singularities and then the geometry 
of the curved bucket is determined through the evaluation 
of the “asayaee Anive Sica Gr. 


Consider the hodograph-plane in Fig. 3b where 


pe eee (10) 


the fluid stagnates at A, accelerates along ANB in a manner 


specified by the hodograph-plane, leaves the bucket at B at 


cal 








fa) z-plane (b) ¢-plane 





(Caro Pp hame 


© 1/k = -k 0 : 
oo 
Be c D CA A 
(d) t-plane 


Figure 3 Transformation planes for the modified-hodograph method 
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moore Of 7, and eventually acquires an angle of deflection 
B. Thus, the outer unit circle in the hodograph-plane repre- 
sents the free streamlines of constant velocity. As stated 
earlier, the fluid is assumed to be inviscid and the real- 
mora ellecus such as boundary layer, entrainment, etc. are 
memored. 

Clearly, the arc ANB in the hodograph plane may be chosen 
Memrepresent differeyr velocity variations along the bucket 
and different angles of departure at B, e.g., the curve AE. 
mienough the analysis will be relatively more complex, there 
MmemeerO Oe sure, no conceptual difficulties and differences 
between that presented herein and other cases. 

The hodograph-plane may be transformed to an S-plane 


(see Fig. 3c) through 


Sees 
2 a1. (11) 


wMmasune interior of the polygon in the S-plane may be trans- 


formed into the lower half of the t-plane by 
el 
S == int ie) 
T 
which is obtained through the use of the Schwarz-Christoffel 
transformation. 


The flow in the t-plane is comprised of two sinks and a 


pemmeee. tence the complex potential is given by 


29 





NEXG 


Whe —l_fin(t+k) in in(t+z2) So @ Mbec ee fa 3) 
Writing 
= eee 
az = ae r aw 


morobtvain the geometry of the bucket, one has 


= Gh Cin : 1 ee 
a nel Yn t Cee 2 ae eal 
k 


] at Cl 


ever Separation of the real and imaginary part, one obtains 


Cums) 
> eae k 
q = = In <a o Cis) 
(t+1) 1 
end 
oe 1 2 
pee! ee eq) ea? Oe (16) 
W aie 


The deflection angle 8 is related to the parameter k through 


SCE (ape) G® 
1 -cosB ; 
which is obtained by writing z = ee iisercing Le in 


Eq. (11) and combining with Eq. (12). 
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iment oregoing completes the formal analysis of the problem. 
ime Ducket f@eometry may now be calculated for a given deflec- 
mom angele 8 through the use of Eqs. (15), (16) and (17). 
Mie results obtained through a straightforward integration 
scheme are shown in Fig. 4 for a set of representative jet 
G@erlection angles. The coordinates of the 50 degree case 
are shown in Appendix A, and the computer program in Appendix 
EeeeoulLiice it to say that the thrust acting on the bucket 
may easily be determined through the use of the momentum 


equation in terms of the deflection angle B. 


Pee GEVI-CIVITA'S METHOD 
This method, which is fully described in several reference 
texts (e.g. Birkhoff and Zarantonello [14], Milne-Thomson 


(4), consists of the definition of a complex function w as 


eae tect) 
J 

where 6 and g are respectively the direction and the magnitude 
oeume velocity Wee cose. mor wee eons tant vVelociuy along 
the free streamlines. The w-function is expressed in terms 
Sieempolynomial in such a manner that w is real over the free 
Streamlines and complex over the curved obstacle. The shape 
Sime Curvature of the obstacle is defined in terms of the 
PCetmmrcients of a polynomial representing w. 

Even though this method was originally devised to deal 


with curved boundaries, as discussed in Section I.B-5, its 


Be 





PmeeecaclLOon has been mostly restricted to boundaries with 
straight-line segments (Cisotti [15]), because of the ana- 
maereamecti ry culties encountered in its application to curved 
Mewmdaries. Wu [16] has employed it in the study of flow 
eeout fully—-cavitating, slightly-curved, hydrofoils. 

The physical z-plane (see Fig. 5a) consists of a jet of 
width ed and velocity _ impinging upon a curved, symmetric, 
two-dimensional bucket BAB’. The complcox-potential function 


W is given as usual by (see Fig. 5b) 
W= 6 + iy (15) 


A straightforward application of the Schwarz-—Christoffel 
transformation establishes a relation between the t- and 
W-planes. Thus, one has 


V.d 
C= as Pinecieatjer in (+t) | Galop 


Because of the symmetry , t is taken equal to tl and -l at 
C and C! respectively and zero at A. The parameter k is 
assigned to B and, as will be shown later, it is uniquely 
Bowaved to the deflection angle 6. 

The upper half of the t-plane is transformed into the 
inner region of a semicircle of unit radius (Fig. 5c) through 


the transformation given by Levi-Civita [12] as 
4 ts i 


Se 








fa} z-plane (b) W-plane 
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(c) t-plane (d) {-plane 


Figure 5 Transformation planes for Levi-Civita's method 
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The circular are BAB' in Fig. 5c represents the solid boun- 
Sime eover Which the direction of the velocity vector is 
assumed to be prescribed and diametric line BCDC'B' repre- 
sents the free streamlines where the magnitude of the velocity 
is known. 

timemrelavionship between the points in the z, t and 


eplanes are given by: 


ardl 


Bomnu 8°, EC = -K; ct 


peimt A , t= 0; CF i 


(21) 
Doin Dt +o, c= 0 


Ommuas: yee io = kK , C= -l 
The W-plane is related to the z-plane, as before, by 


~ v an = ye G22) 
Where q and 6 represent respectively the magnitude and direc- 
Ptommot the velocity vector. 

One can now define an w-function which is quite similar 
uO tnat given by Planck as 


1 QW) . ae 
7 an? Cte 2 (F ) C289, 


J J 


wW = i in (- 


Now solving for z in the above equation and expressing all 


the variables as functions of Zf one has 


a2 





Z._ 1 -aiwls) dw 
= Bare JG a dz (24) 
J 
Defining 
CF r eld 


pe@esioving that over the bucket only the argument and along 
the free streamlines only the modulus varies, the combination 


preeecis. (19) and (24) yields, 


. o : 
z e » . fertwlo) 2k ooae a) ae (25) 
1 -—- k” cos a 
and 
Ke, 2 Ln 
ee fu(p)cor 7 re? 
a a [ } ar (26) 


2 
<(e + =) 6 - 


Equation (25) gives the coordinates of the bucket and Eq. (26) 
the coordinates of the free streamlines. 
The function w is assumed, according to Levi-Civita, to 


be of the form 


aoe 
2 er 





—— 2 3 
w(t) = i ln AB Sar Gk (6 As Ae me a eee C27) 


2° 


where the logarithmic term can be identified as that corres-— 


ponding to a free jet impinging normally upon a flat plate 





Mersoutrili7)). The additional terms in the polynomial whose 
Seerlicyenves are to be prescribed later modifies the expres- 
rom tor ow and helps to curve the solid boundary. 

Because of the symmetry of the flow in the physical 
mame, the coefficients Aos Ao. Ays +++- are not AGO SOE. 

Noting that over the solid boundary r = 1 and over the 
free streamlines o = 0 oro = wt, Eq. (27) may be reduced 
over the bucket to 


io 130 


te aoe hae ee 


We] 2 im Leos Oo) 


3b ST) 
Ai +f Sea no _— ee e 8 


hola 


5 
(28) 


where plus sign is to be used for tT > o > w/2, and minus 
sign for 1/2 > o > 0. Over the free streamlines, it reduces 
iGO 


i“) = — tan ~ a arta pe 


tar... (29) 
The coefficients ay> 23s Ap» ieee May NOW De chosen tO. obtain 
different bucket shapes. In the present study, w is chosen 


to Sataisity the following conditions 


w(+1l) = =n 
misty S =o alec 
(30) 
w(-1) = 4 
wlo) = 0 


att 
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Meese conditions simply state that the jet departure angle 
femme: Jip Of the bucket be equal to tn, and that the bucket 
be continuous in curvature at the stagnation point A and 
Symmetric with respect to the jet axis. 

There are several ways to select the coefficients as: 
In this study a Fourier-series representation was used 
Peeemari ly because of the versatility afforded by it in 
Gepresenting different bucket shapes. 


Equation (28) can be written over the circular arc AB' 


as 
wo) = if{ln Peenses ay Sino + a2 SWigee Sieh eec 25 Syiliel oyekneen |) 
(se) 
TI 
+{- 5 + a, cos o + a2 COs S008 as COS Ole 


Mgeereal part of Eq. (31) represents the direction 6 of the 
Tmeeciuy Vector over the bucket. Thus, it is obvious that, 
one can prescribe the curvature of the bucket or the modulus 
ei tune velocity eter by properly choosing the coefficients 
ay through 2. 

Let @ be represented by the lines shown in Fig. 6. Evi- 
dently, by choosing sufficient number of terms, one can re- 
present the variation of 6 with o shown in Fig. 6. In doing 
so it is important to make ® symmetric with respect to the 


6 axis since the series are expressed in cosine terms only. 


Rewriting Eq. (31) one has over the arc AB! 
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Figure 6 Approximate representation of O(c) 


3 


ya 





COS @¢ 


ac) = tn ic. ¢ h(o) + 6(a) (32) 
where 
N 
Mg) = Y @. sin no (33) 
a n 
n=l 
and 
, ON 
6(a) = Bee ee COs" MG (34) 
= ri 
n=l 
where 
Ae es See atthe tae ae [eos no. - cos no. |] (35) 
n 2 al 2 
(o4 -- o5)n 


aa Which ony and 5 represent the two parameters shown in 
Fig. 6. The evaluation of the Fourier coefficients a alts! 
Emeowd in Appendix C. To complete the analysis one needs to 


arte) 


evaluate the terme Wile ampears: in haq..9(25) in terms 


er Ba. (32). Writing 


eiulc) Meh a Ae (cy) ar8(o) 
cos Go 


and separating the real and imaginary parts, one has 


etw(o) . L+sing ehlo) (eo; g(g) + 4 sin 6(o)) (36) 
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few combining Eqs. (25), (33), (35) and (36), one finally 


aocains 
x a g ~h(o) 24° sin of] + sin og) 
ef we Sin 00) Saniora amr eae ss 0 
tive l1- k eos a 
(37) 
and 
O Ces 
d : 17 2 1- 1 ae 
ip) 


Equations (37) and (38) yield the coordinates of each 
meemens Of the curved boundary through the use of the appro- 


Mmatace values of go, ando.. This in essence completes the 


Ht 2 
application of the Levi-Civita's method to the analysis of 
tae jet deflection from curved obstacles. In the following 
Mae amplications of the variations of 04 and 5 will be 
@ascussed. 

The variation of 6 with o shown in Fig. 6 may be written 


for each interval as 


Q(o) = =n OF-moe< ort 
el 6) Tt el 
6(0) = 5 ( ) - By (2 + ) O4 <0 < 05 
rey tee al oA 
Q(o) = - 3 Sn <0 < z (39) 
6(o) = - 5 5 <o< 17 - 05 
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imeewill be shown that this particular functional depen- 
memce produces bucket geometries of the form shown in Fig. 7. 
Rewriting Eqs. (37) and (38) for the particular example 


under consideration one has 


Xx 
+s = © (40a) 


and 


Uae S if a sy, Ke See aches oO. 
TT ras ee ree 


1/2 ey ke tres ae 


dos (40b) 


mor the interval of 7 > 5 205; i.e. for the segment AP, 


(see Fig. 7) 


O Eis os o- 3G 
x Z é f Zot (o) ius os i Ae + ae o) cos(% —— : pee 
05 1 - k' cos a 2 A 
Cale) 
and 
0 2 oe. Oo - oO 
M = 2 r en tho es Eee cee OY (3 —— : lec Y, 
D5 1=—- k cos oa 2 uh 
(41b) 
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Menmeevne interval 0, > 6 > oO i.e. for the curved segment 


2 Hee 
PQ, and finally 


0 ies 
: _ é rp ents) k” sin ss aE sto 2a ae (42a) 
O5 Dea e eeksl ire) 4 
and 
Y = 
a G (42b) 


for the interval o, > o > 0, i.e. for the segment QB. 

The coordinates of the free streamlines may be obtained 
by integrating Eq. (26) through the use of Eq. (29) and the 
coefficients a. Mivemmeby vid. (3/)en Pineally., the det lection 


angle is evaluated by replacing in Eq. (29) the corresponding 
maluve of r at the point C, i.e. r= = (1 - v1 - Ke yee 

The foregoing equations have been numerically integrated 
(the computer program is given in Appendix D), and the re- 
Sulting family of bucket shapes are shown in Fig. 8a. A 


= 0.31416, o, = 0.47124, and k = 0.946 


special case for o 


ale 
fioeshown in Wig. Ob. 


C. RIEMANN-HILBERT METHOD 
The method consists of the transformation of the physical 


plane into the W-plane, where W is given by 


W= ot iy 
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Pmemune tCransiOrmation of the W-plane into the upper half 

of a t-plane through the use of the Schwarz—-Christoffel 
transformation. Then the flow in the t-plane, which is 
composed of segments of solid boundaries and free stream- 
lines are related to the Riemann-Hilbert transformation as a 
well-posed mixed boundary value problem. Then a certain 
Shape function f(t) is ‘assumed and the characteristics of 
the flow in the physical plane are evaluated through 
micesracvtion. 

Mie cevaills of the technique will become clear as it 
mmeopplied cto the determination of the deflection angle of 
a two-dimensional jet impinging symmetrically upon a 
two-dimensional bucket. 

Figure 9a shows a jet of width 2d impinging upon a curved 
bucket with a 180 degree angle of departure at the upper and 
lower lip of the bucket. 

It is apparent from the application of Bernoulli's equa-~ 
tion that the velocity is constant along the Tree ssuriraces: 
The complex—function W, for the upper half of the physical 
plane, is represented by an infinite slit as shown in Fig. 9b. 
The W-plane can in turn be transformed into the upper half 
of a t-plane through the use of the Schwarz-Christoffel 
Vranis sora tole sine, surarenulorward application «of this 


transformation yields 


A Aero! 
wit) = “~~ int (42) 


TT 
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where we is the velocity along the free surfaces. In 
accordance with the rule of this transformation, three of 
the t values are assigned arbitrary values (t = 1 at A, 
t = to at C and t = 0 at D), and an unknown parameter k is 
assigned to the point B. 

The complex-potential-function W is related to the 


physical z-plane by 


= oo (43) 


fmere q and @ represent respectively the magnitude and the 
Oivecs.on Of the velocity vector. 

it is convenient at this point to separate the magnitude 
meemcne Glrection Of the velocity through the application 
om Planck's logarithmic transformation given by 


aie O (44) 


WwW = In on 
J 


Pomoinines Eqs. (42) and (44) and solving for z in terms of 


t, one has 


a 1 -w(t) dat 
= a Be o (45) 


O(N 


It is apparent from the foregoing development that the real 
part of the complex-function w(t) is zero along the free 


streamlines where q = Vee Thaws ; 


48 





Refw(t)] = 0 for ‘eee and arcane 


The imaginary part of w(t) which represents the direction 
meethe velocity vector is zero along the line AD where 
Wea t < |. Over the bucket, however, both the magnitude 
and the direction of the velocity are variable and are an 
meaenown function of t. Let this function be denoted by f(t) 
mmecne interval 1 < t < k. 

The foregoing conditions define a well posed Riemann- 
Hilbert mixed-boundary value problem where the real and 
imaginary parts of an analytic function w(t) are alternately 


rad@wr alone the real axis, 1.é., 


Re(w) = 0 (eae meille 

Im(w) = -f(t) ee pies kc (46) 
Im(w) = 0 Gece 
Re(w) = 0 t < 0 


ewe cdc owecolurion of the Riemann-=Hilpert problem is well 
Mmioune«chenesandanorl. . 10) and Larock [9 ]). 


a Dein pre mone” Wirites 


Wiciea) ullsl Ceo aien (47) 
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where H(t) is the solution of the homogeneous—Riemann prob- 
mame c. for £(t) = 0) as it satisfies all the other condi- 


tions listed in Eq. (46) and may be chosen as 


H(t) = y -t(t ~ k) (48) 


in passing it should be noted that the particular form of 
mec) which constitutes a solution for f(t) = 0 is not of 
special importance as shown by Song [19] and that uniqueness 
Sethe solution does not depend on the particular value of 
met). 


The function Q(t) may be written as (Larock [9]) 


+00 
at) == f  SOM! an (49) 


= OO 


Wembining Eas. (46), (47) and (48), one has 


a(t) = Yet =) __ =f) dn 
and td tn = 


(50) 


The remainder of the analysis depends on the selection 
of f(t). As pointed out earlier, f(t) would have been known 
had the shape of the solid boundary been known a priori. The 
present analysis, however, is an indirect approach and seeks 
the shape of the boundary in terms of the boundary conditions 
assigned and the deflection angle assumed. Thus, it is 


apparent that one may select f(t) as an arbitrary function 


DO 





Mmmee provided that the boundary conditions at the stagnation 
meepeeeand at the lip of the bucket are satisfied, and obtain 
meramily of bucket shapes. In the following, f(t) will be 
taken as a linear function of t partly for the sake of sim- 
meerei Ly and partly in order not to obscure the AP cers 
of the method by introducing rather complex and unmanageable 


mavegrals. Thus, writing, 
f(t) = a + bt es) 


ene choosing a and b in such a manner that the slope of the 


bucket be + 1/2 at the point A and + m at the lip, one has 








f, nlg eee 
Ca = 1 
and 
gene 
a ea 


Pmoercing Ba. (51) in Fa. (50) one obtains 


_ -(a a bn) dn 


lig Syma etie 


—~t(t — k) 
Tt 


w(t) = (52) 


The integration of Eq. (52), the details of which are 


presented in Appendix E, yields, 


le 





fOtemUresmne calla, b> ik 


and 
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aati le hens Sea 


Citbe) 
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t(k-t) cos + | 2 | 


w(t) = k 
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WimesdetTiection angle 6 can be found by letting t¢7* -~ in 
Eq. (53), or by letting t > » in Eq. (50) and then performing 
the integration as shown in Appendix F. The deflection angle 


8 is found to be uniquely related to the parameter k as 


B= tt oa | ( - 20087? 2K + yx -2| (55) 


To obtain the shape of the bucket and.the actual location of 
the free streamlines, Eqs. (53) and (54) are replaced in 
Eq. (45) and integrated. 

The integration must be performed numerically using the 
method described by Hartree [20] because of the singularity 
at the point = Jeo hans. Den pointed out thatthe 


Singularity at t = 1 is an integrable singularity of the form 


py 





po ox 
0 Waa 
Memacvcvempt iS made here to perform the numerical integration 
of Eq. (45) even though the procedure is well known as shown 
in one application by Tinney [2], because the method previously 
described, namely, Levi-Civita's method, is very similar to 
the one presented herein and has the advantage of being very 
flexible to obtain different bucket shapes and presents no 
eamezularities for the integration. 
eee CRITIQUE OF THE ANALYTICAL METHODS FOR 

TWO-DIMENSIONAL FLOWS 

The past as well as the present analyses have emphatically 
eiietouraved clao almost all of the difficulties of the 
analysis of jet deflection from curved obstacles are ascribable 
to the determination of the initially unknown free streamlines. 
irevtew Of Gars fact a direct analysis of the problem, that is, 
wae Cetermination of the jet deflection and the position of 
the free streamlines via either purely analytical or partly 
analytical and partly numerical techniques is at best an 
exceedingly complex proposition. As propounded by Birckhoff 
et al. [14] and other hydrodynamicists what one can hope for 
ToeoMeniCireChrsOLUer On OrmLie Droblem;othal 1s, obvalning 2 
family of obstacle shapes which will yield the prescribed jet 
Gelleevion angle and provide a familiarity between. the apenee 
SO Obtained and the range and change of direction of the free 


parameters involved in the analysis. What has been said so 


Be 





far holds true for both the plane and axisymmetric flows. 
Should one adopt a purely numerical procedure (finite differ- 
ence, relaxation, finite element, Spe) ate sw Unem. and sony 
men that it is possible to begin with a prescribed obstacle 
shape and to approach a unique solution through successive 
meeratvions. 

Miewhat follows, attention will be focused only upon the 
evaluation of the individual as well as relative merits of 
the analytical methods dealing with the deflection of plane 
jets. 

The hodograph-method developed herein is capable, within 
mem iimitations of the potential-flow theory, of analyzing 
the flow in curved obstacles and generating a family of 
practically suitable bucket shapes. Furthermore, this method 
Meeethe inherent capability of being extended to the genera— 
tion of even a larger family of obstacle-shapes through the 
use of distributed sources and sinks in the hodograph-plane. 
mem tnough this concept has not been pursued further herein, 
there are no conceptual Gi¢trculties ine doing so. 

The method devised by Levi-Civita and successfully em- 
ployed neve requires the judicious selection of the varia- 
@lton of 0 with o and the use of a suitable series to enable 
one to perform the integrals leading to the coordinates of 
the obstacle as well as of the free surfaces. The deflection 
angle may be determined either by calculating the asymptotic 
slope of the free streamlines or by finding the pressure 


distribution over the bucket and using the equation of 





Mometnvum. The emerging fact is that a brief familiarity with 
torreemethod enables one not only to generate a family of 
meee le curved obstacles but also to round the corners of 
buckets composed of straight segments. For example, by 
letting (oO, - o.) + 0 one finds that the bucket shape 
Pemeraved by this method exactly approaches that previously 
peeved by Sarpkaya [1] (see Figure. 8a). For slightly 
larger values of (o, - o.,) the corner ®@20 Smee Fig. 7) is 
Moree d, thus enabling one to evaluate the effect of the 
deviation of the segment PQ from a sharp corner on the 
merece vion angle. in fact, the results show that for a 
given b/a and a/d ratios (see Fig. 1) the effect of the 
curvature is relatively small as anticipated by Sarpkaya [1] 
Gaeenme erounds that the flow of a real fluid near the sharp 
@ormers form nearly stagnant regions comprised of small 
Corner vortices and provides a slightly rounded streamline. 

Finally, it should be noted that Levi-Civita's method 
has been rendered practicable through the introduction of 
the dotneuateneermae imuvo ChevanalySis.0) Lt may be ECasily 
extended to the analysis of the cases where the lip angle of 
enesebUCKeL 15) 0uher than 1 or the bucket has a cusp on its 
axis of symmetry or to the cases where the impinging jet is 
not free and emerges from a prescribed Nozzles (Sveh shapes 
are known to be used in the design of spillways, Pelton 
wheels, two-dimensional cascading thrust-reversers, etc. 

As to the Riemann-Hilbert method, it may have already 


been evident that it is in essence a mathematically more 


oye, 


re 





elegant and sophisticated version of Levi-Civita's method in 
meoscriDping Certain functions. I1t is for this reason that 
no particular solutions have been presented and only a linear 
variation of f(t) was studied. Needless to say, a quadratic 
or cubic variation of f(t) could have been employed. 
Evidently the resulting integrals would not have been 
miverrable in closed form. Finally, it should be pointed 
euv that the RiemamiMilbert method is perhaps the only 
method aside from purely numerical methods with which one 
PmmrncOrpOrave into the analysis the efiect of body forces 
meen as Sravity. 

ineelosingerynie discussion Of the methods dealing with. the 
G@etlection of plane jets from curved obstacles, one must once 
more emphasize that all of the foregoing analyses and solu- 
tions deal with plane Laplace-field problems where the real- 
fluid effects must be out of necessity ignored. Because of 
The limitations imposed upon man's understanding of turbulence, 
the development of laminar and/or turbulent boundary layers 
under a nonzero beeenee gradient along curved walls, entrain- 
ment between the get and the surrounding fluid which may or 
may not have the same physical properties, and finally the 
Coanda effect (the tendency of jets to attach to adjacent 
surfaces because of the entrainment deprivation), the real- 
fluid effects on jet deflection cannot as yet be analyzed. 
It is because of these reasons that the analysis of the 
PEeopleciaWwiErine Lae aseumpeLon Of an inviscid 1luid becomes 


more important for it is only through the comparison of such 
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solutions with the experimantal results that one can assess 
mm@eeaesree of significance of the real fluid effects on the 
Beverning parameters. As stated in the introduction and as 
verified by numerous experiments (Robertson [13], p. 626), 
meee phenomenon of jet deflection is governed primarily by 
mmessure and inertia forces and that the reai-fluid effects 
meewei Minor significance in determining the forces acting 
emetene body and the deflection of the jet.. These facts will 
memouressed again in connection with the discussion of the 
em@auyctoical and experimental results for axisymmetric jet 


eoerection to be subsequently presented. 
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Pi ANALYSIS OF AXISYMMETRIC JET IMPINGEMENT 


A. INTRODUCTION 

As stated in the general introduction, the two- 
Gimensional counterpart of the jet-deflection problem has 
been treated by several investigators through the use of 
numerical and approximate analytical methods. Additional 
techniques based on the modified-hodograph method, Riemann- 
Hilbert integral transforms, and the Levi-Civita's method 
have been successfully developed and employed in the 
mreceding section in the analysis of the jet impingement 
upon two-dimensional curved buckets. 

The three-dimensional counterpart of the jet-deflection 
problem has not yet been solved in any generality. Attempts 
feom1ormulate an exact solution have been mostly unsuccessiul 
evem Lor the axisymmetric inviscid flows with no body forces. 

The case of a circular jet striking a plate normally was 
analyzed by Schaach [21] using approximate methods similar 
to those of Trefftz [22] with successive adjustment of the 
free streamlines and through the use of the Fredholm integral 
equations. These results are in good agreement with those 
obtained experimentally and to the surprise of the many 
investigators, they are almost identical to those of the 
two-dimensional cases. 

Other noteworthy contributions to the analysis of the 
Jev-elfiux from nozzles and orifices with straight boundaries 


were made by Southwell and Vaisey [23], Rouse and Abul-Fetouh 
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mer), Garabedian [25], and Hunt [26] through the use of the 
relaxation and finite difference techniques, and singular- 
integral equations. Schnurr et al. [27] used the relaxation 
method to analyze the turning of two-dimensional and axisym- 
Mmewric jets from curved surfaces where there was only one 
free stream surface; i.e. the jet was assumed to leave the 
deflector exactly parallel to the tangent at the lip of the 
Semmector surface witb .aeuniform velocity distribution and 
the consequences of the difference between the actual deflec- 
tion angle and the said tangent to the deflector surface was 
taken as a measure of the "spillage" and expressed in terms 
of a "turning-effectiveness coefficient" determined experi- 
Mearcally. Thus, this analysis does not constitute a solution 
to the problem under consideration. 

Jeppson [28] proposed an inverse formulation of the 
Mme@oten for the flow exitirane from a circular orifice by using 
x and r as dependent variables and ¢@ and Ww as the independent 
variables. Through the use of this technique Jeppson [29] 
also analyzed the flow of a jet from a convergent nozzle with 
emeved boundaries. 

All of the investigations described so far have attempted 
COupPOvIGembercer SOlULLONS vO the various econfigurations “of 
the orifice flow. ~These techniques suffer in general from 
convergence and accuracy problems and nearly all resort to 
simple trial-and-error procedures to locate the free surfaces 
and to satisfy Une= DOUNGary CONCItiOon that the ‘tree ‘surtaces 


be streamlines of constant velocity. Suffice it to say that 


a9 





meme vnod was needed which could yleld solutions of a 
Seeescribed accuracy for a wide variety of fluid flow problems 
imeolving Dirichlet, Neumann, and mixed boundary conditions. 

Zienkiewicz and Cheung [30] proposed in 1965 the applica- 
men OL the finite element method to the solution of field 
problems involving the equations of Laplace and Poisson. 
pamce then a significant number of applications of the finite 
element method to fluid dynamics has annveared in the litera- 
ture (see e.g. Norrie and de Vries [31] and the references 
cited therein). 

The method has been recently applied to several jet-efflux 
problems involving only one freestream surface and relatively 
prem jet contraction by Chan and Larock [32]. Suffice it to 
say that the finite element method has proven its versatility 
and applicability to a wide variety of two-dimensional and 
axisymmetric flow-problems. The truly three-dimensional 
peep lem Such as the oblique impact of a circular jet on a 
peeane SUriace, With or without gravitational forces, just to 
name one, offers considerable challenges and demands additional 
Sepatsculcavion. 

The special Case sol “Une Gravivationalvioree acting along 
the axis of symmetry may easily be analyzed by modifying the 
RPevUdaary concieutoncesuch thas the velocity alone the tree 
SUrrace  Variesy according COsthe Bernoulli theorem. 

This part of the investigation encouraged by the advent 
and the continued success of the finite element method is 


devoted to a determination of the angle of deflection, the 
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meere1on OL the freestream surfaces and the velocity and 
Pieesoure distributions caused by a finite curved boundary 
memeeneral and a hemispherical boundary in particular placed 
eeometrically with respect to the axis of an axisymmetric, 
fieeesocia jet issuing from a nozzle. This problem has not 
been previously analyzed through either exact or approximate 
methods. 

As stated earlier, the deflectionwf a jet by a solid 
moundary 1S well suited to potentiali-flow analysis because 
of the dominance of inertia and pressure forces in the 
establishment of the flow pattern. The result in general 
mepresents asymptotic values which are approached as the 
effect of secondary variables such as entrainment, boundary 
Mayer, compressibility, jet attachment to adjacent surfaces 
(Coanda effect), etc. are decreased. Specifically, if the 
total angie through which the jet is deflected is determined 
memecoealizead conditions, then the principle of impulse and 
ieee uni Callebe Used FO compute pressures , rorces, or other 
Gynamic characteristics of the system such as the reverse- 
Mimchiotrraulo. Mins tmis Manner one can nov. only optimize the 
geometry of the system for optimum performance, measured in 
berilns. Ol reversea thrust, but also can assess the signiticance 


of the real-fluid effects, cited above, on its performance. 


Bes ANALYSIS 
The problem may be defined in general with reference to 


Figure 10. An axisymmetric jet issues from a nozzle BC 


Gal 








Bieure 10 “General case of axisymmetric bucket an nozzle 


with an entrance velocity V, at AB. The jet impinges symmet- 
rically upon the curved bucket LKH and then the deflected jet 
eaves av an angle 6. The upper and lower surfaces of the 
Hew. namely. CDE and HG, asyvinouoewcaliy converge: Corinne 
eee ope the jet in accordance with the conservation of 
mass. 

Lembo apparentmurom Mie.) tO thay <i or ia elven nozzle and 
bucket share, the deflection angle varies with the diameters 
of the bucket and the nozzle and the separation distance 
between the Tips ol vine nozzle and the bucket. In other 


words, 8 can be expressed as 
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8 = wae 3 S/S) C56) 


mms, the determination of this functional relationship 
Sematvitutes the basis of the problem. 

waemiluid as assumed to be inviscid, incompressible, and 
free from body forces. Consequently, the problem may be 
fommulated in terms of either the velocity potential $6 or 
fmgemovookes stream function W. Here the potential function 
memenosen as the unknown field parameter. 


The governing equation is then given by 


iameidcl X and Pare the axial and radial coordinates. 
Because of the axial symmetry there is no variation in the 
BeamarOnal direction. This equation is to be solved with 
The boundary conditions which simply states that the normal 
Component of velocity be zero along both the solid and free 
Buia Ces. - i.e. tHe bree. surtaces Dev slreamlines OL -Constany 
velocity. 

HOr yan axisymmeuric flow, the solution to the Laplace- 
field equation satisfying the specified normal-—-velocity 
boundary condition (3¢/3n)* is eiven by the admissible 
function $ which minimizes the functional (see e.g. Zienkiewicz 


and Cheung [30]) 


7 OO dd ,2 od,a 
I(¢) = ot . ($2) + (<) jJrdxdr - 207 ; o(S2) rds (58) 
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maewietch A is the half of a meridional section of the flow 
meoerc 1S a2 portion of the curve bounding this area where the 
mormal derivatives are prescribed. The first integral in 
Eq. (58) represents the kinetic energy of the fluid within 
mae entire control volume and the second integral represents 
twice the work done by the impulsive pressure pd on the 
boundaries in starting the motion from rest (Milne-Thompson [4]). 

It is apparent from the foregoing that in the finite- 
element method the boundary conditions are an integral part 
eaethne functional I(¢) and will not have to be dealt with 
separately as in the case of finite-difference or relaxation 
techniques. This is, in fact, where the beauty and the power 
Ot the finite element technique lie. 

The approximate minimization of the functional I(¢) is 
eecomplished as follows. Firstly, the field of interest is 
divided into M triangular elements and a functional I(¢) is 
written as a sum of M element-—-functionals 1° (6) each of “whlch 


Hecethe form of Eq. (56), i.e. 
“re 
CG j= ee ee) (59) 
il 


Figure 11 shows the numbering system of the corners and mid- 
points and the corresponding area coordinates GC. = A. /A” 
in which A, is the area of a subtriangle and Af& is the area 


OEeenexcHyuurcotrianeile, Sinenmit 1s clear that aT + 5 + 53 ole 
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Figure |! A triangular element with area coordinates 


peconcly, a polynomial ot sufficient order is selected to 
represent the variation of the potential function © within 
Cae ec remenu. 

Apparently, a linear Variation of © leads to constant 
velocities along a given line within an element and to 
unacceptably large errors in the velocities between neighboring 
elements. <A second-order polynomial, however, allows at least 
a linear variation of velocity components. 

The use of higher order polynomials may increase the 
accuracy of the results and the convergence of the minimization, 
BUG Aly co also lead to extremely tedious arithmetical manipu- 


lations. For these reasons and following Felippa [33] and 
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Chan and Larock [32], a second-order polynomial has been 


@mosen. in terms of the area coordinates, this leads to 


ae ae en Cn Pe) (60) 


where b5 eee UMe COC Iclemes Ol the polynomial. and corresponds 
moeecne Values of aS Hee cMemlOGesn .¢ “10 ls 2NpOrtane vOunoue 

whe t b5 Wteu De srerarded as constants during the integration 

of Eq. (58) in each element. 


The term X, may be written after Felippa [33] as 
X_ = [eo (20y-1),050205-1) 0 3(20 371) Meta sNoob354t S32) (61) 


Mae velocity components u and v are given by 


n= 2 
OX 
and 
- 9 
Mf or 


Men the velocity in each element can be written as 


| dX 
CL. al ate 
bas os a ‘et i everioe (62) 
ler 
e _ .e af are 
V = Ps Oa ea = ] to 6) (63) 


The evaluation of Eas. (62) and (63) require the derivatives 
Of ye with respect to x and r or the evaluation of terms 


such as 
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me oc. =) oT b 
OX OX CA 
mec —1) ava b 
mee’ °2 lee (1 poe oo hee 
SS = CU - 1) 3 = (405-1) — = T 
yx o OX O one e 
ee(2ec.-1) aT b 
a3 3d = (ue 1) —3 = (4z,-1) 4 = 18 
OX OX OA 
a(4e.c0.) aT of 
es yr ae = re 
ax ls Sears ma ere Dae ()b, + S,b,) = T, 
0(4e.7.) ara ora 
er eb) Ds oe = me 
OX 7 oo OX + AC. OX ape (5b. 4 Gb.) ae 
Oe HG it) of oT 
aaa 225) i ee eae = me 
ax Reena Eee) OK ae (5,6, + o,b,) = T¢ 
in which 
JE Ease 
ee) a 
and 
D. = rs - i 


The foregoing closely follows the notations developed by 


Meslbatjo)oysa less ial 
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7 


ieee evideny thavethe derivatives of x with respect to 
mea Casily be evaluated by changing x to r or bd, tO a. 
given by 


pmeecenoting the resulting array of qT. by qT, - Thus -Cne 
can write 


(4c, - 1)b. 


jt a a eee (no sum on i) (65a) 
A ZN 


eiGbeee at b5o4) 





e z niga 
eee one (65b) 
and likewise 
ie gan 73 l)a, 
TT. = Cacesulneone 1) (66a) 
ae e 
OA 
pe SCR ew aan) 
S42 2 4 (66b) 
it+3 Bie 
In which 
a. : aan 
one ar 
ana 
r= EST a P5o5 a 3 P23 
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iiommaninazavton- ot the functional 1(¢) in accordance with 
the Ritz method requires that the partial derivatives of 


I(¢) with respect to >. be zero. In other words 


Dee 
at {8) = 9 | (67) 
o>. 
which may be evaluated in a manner quite similar to that 


Misea in Eq. (64). Thus, Eq. (67) reduces to 


al" (o) _ ee. ee, ,ée ao\a - 
242 = 2pm 7 (TET, +T5T5) os rdA - 279 : x, (Sh) "ras = 0 (68) 


or CO 


pie Col) ois) eae e 
ee Ki 5% =/B l= 0 (69) 


265 

in which Ky and By represent, respectively, in the structural- 
mechanics terminology, the element stiffness matrix and the 
load matrix for a triangular element. In fluid mechanics 
terminology, they simply represent the geometrical configura- 
jcnkol@enond the erid within the domain and the boundary conditions 
imposed on it. These coefficients have been previously 
obtained by Chan and Larock [32] and are tabulated in Appendix 
G. 

In fluid-flow problems where one or more parts of the 
boundary are to be determined as part of the numerical 


analysis, the solution of the system of equations represented 
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by Eq. (69) through the use of the Gaussian elimination 
procedure requires the use of a suitable (i.e. accurate and 
convergent) iteration scheme. This important phase of the 


analysis will be described in the following. 


C. ITERATION SCHEME 
The iteration scheme consists of the designation ofa 
meving and fixed grid, the adoption of a scheme for assuming 
the initial position of the deflected jet, and the successive 
Scbsemene Ol ie Iree suriaces until. the boundary conditions 
ame satisfied. 
im, elesw eau Om.or la Grid 

Evidently, the minimization process cannot be carried 
MUEeWiol a fixed Grid as it would be in a structural system 
subjected to small deformations, primarily because the 
ime rdieMOoote rom Ol. Lhe tree surtace 1S: not <a priori, known, 
and the initially-assumed boundaries may be too far off from 
iicim tina he DOs tOns. imi view rot the, foregoing. the Tiuird 
domain was divided into a moving and a fixed grid (see Fig. 
12). Rien ierxccmorid: Occuples that regvon-O1 “the flow which 
MelOu.eligeeme tinal analysis.) to be inversecvted by Une inec 
Surtaces. sIhis does not impose any restrictions on the anal-— 
ysis and requires no more than a passing familiarity with 
wie eleMmehcaryveprmiciples Of jen deflection from Solid 
boundarice.. o1oulcdepme DOUuNn@dary Ol ‘the fixed ‘rid be anver- 
sected by the free surface during the iteration process, a 


new fixed grid has to be assigned. 
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The moving grid is comprised of triangular elements 
whose nodal points may be moved along prescribed lines prior 
mowtne cCOmMmencement of the next iteration. This procedure may 
meee be described with reference to Fig. 13. At the end of 
emerven iteration, the nodal points A, B, C are moved inward 
or outward, in a manner and magnitude to be described later 
meeemcunie pOSitions of the nodal points F, G, E, H, D, and K 


mee recalculated such that 
jE SEU (Ole alls a eG 
and 


AG = 0.50 AN 


and similarly for BE, BH, etc. Evidently, the coefficients 
Oris ana 0.50 are arbitrary and different values could have 
meen Used provided Chat. the resulting elements maintain a 
miape more or Tess compavible To their original shape and 
iat Tie Tarlo on thie=sildes of a @iven element does not 
Suen tl ecantly nener ee) ie Veg ak ec 

The points intermediate to AB and BC, i.e. the points 
1 through = are calculated at each iteration through the use 
of a subroutine called CURVE-FIT. This subroutine is a third 
order polynomial connecting the points A, B and C, and evenly 
spacing the points 1 through 3 between AB and the points 4 
through 6 between BC. It should be immediately pointed out 
that there is no need for an exotic subroutine to perform 


this procedure and Chat any other judiciously-selected 
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Figure 13 Construction of the moving grid 
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femeynOomial could have been employed. The sole purpose of 
the foregoing procedure is to smooth the boundary between 
two successive iterations and to obtain a better representa- 
mom Of the curved streamlines. Furthermore, this procedure 
mmonaves the Occurrence of localized low or high velocities 
Which could have stemmed from a zig-zag positioning of the 
Besnes such as 1, 2, 3 and thereby accelerates the convergence 
process. In the deflected portion of the free jet (see Fig. 
io) tne points such as P and R are moved inward or outward 
pomclctpaced by the iteration conditions and the point Q is 
iecaced in the middle of PR prior to each new iteration. 
2. Assumption of Initial Boundaries 
iiceassuUlpelOn Ol al initial boundary requires. the 
memotderation of the equation of continuity, some familiarity 
With jet deflection problems in fluid mechanics, and possibly 
mae use Of aderlection-angle predictor, if one is available. 
Ordinarilvy,sa.deilection angle 8 and the length of 
the deflected free jet (HG) are assumed (see Fig. 14). Then, 
o? TR ano LRG equation: Of~ GcOntanuity. 


the thickness d of the jet is calcuated, to a first order of 


through the use of Wb ig 


approximation by 


O 
O 





SO Ror (70) 


Subsequently, the points G and E are located such that F is 
imimermeathauve vO DOuUN. Then the point G 1s connected to H 
(rnewupper dian of the tet) by a straight line, and the point 


tas SsmOoOunhly j@anea Fo the Lower portion or the assumed 


TH 











Figure 14 Nomenclature for a hemispherical! bucket and 
straight-nozzle combination 


beer y. tn passing, 1t should be noted that the approximate 
Moe ol cCcne equation of continuity in sections intermediate To 
GH and ED may significantly simplify the assignment of the 
Piicial position on ED.) For this purpose, all one has to do 
is to assume AMM hOrimevelocity at each. section and calculace 
MiecwuiNUeKness Or Une jet as the section under consideration 
through the use of an expression similar to that given by 
Eq. (70). 

The description of the assignment of the fluid domain 
will not be complete without a brief discussion of the length 
Olav Mewmeczile- and tne deflected jet. Clearly, the sections 


AB and EG must be so located that one can for all intents and 
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purposes assume the velocities at the said sections to be 
uniform. Otherwise, the analysis will not be correct and the 
Peaeevoundary will be influenced by the conditions created 

by the proximity of either the section AB and/or EG to the 
Permmtcowe and fi. To this end, preliminary calculations were 
Carried out with nozzle lengths (BC) as large as br. eVgiell ihc 
Wasceiound that 10 could easily be reduced To oe without 
Materially affecting the accuracy of the velocity at AB. In 
fact, the calculated velocities at AB with a nozzle length 
of Sige were accurate within + 0.5%. Calculations were per- 
formed with a nozzle length of 3.8r ,i.e. about twice the 
mUCKen  Glhanece? . 

The deflected-jet length GH was varied between Oe 
and O.72r. depending more or less on the particular angle of 
Getlection. Surlice it to say that no noticeable change was 
@pserved in the calculated position of the jet. This is 
primarily because of the fact that for a jet length of O.6r, 
Mie ratio Of the jeu lenetvh to final jet thickness is about 
S.0 and that oe aOeae to ai thes previous .calculat tons, and 
exact solutions performed for two- and three-dimensional jets, 
ines el euvboLMsultomesiMpuvCure Unitorm velocity profije in 
aedistence Wess than 3d. 

Seve Vou lon Oust ier mp OuUMCdaInnes 

In the foregoing, all of the necessary elements of 

the analysis and the domain have been provided for the com- 


mencement of the iteration process. We will now describe 
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mew process which eventually leads to the final position 
Oemeecne jet for a given geometrical positioning of a 
hemispherical bucket and a straight nozzle (see Fig. 14). 

To begin with, the free surfaces were assumed to 
be solid boundaries. Thus, the two boundary conditions; 
mamely, constant pressure and zero normal velocity along the 
free surfaces, reduces to a single condition. Specifically, 
the freestream surTeees shoulé@ be streamlines of constant 
velocity. The normal velocity was written zero everywhere 
along all the rigid and free surfaces except along EG (see 
Fig. 14). An arbitrary value of @ = 100 was assigned along 
EG. Needless to say, the final solution does not depend on 
Peers particular value of 9¢. 

To start a computation the velocity V. along AB was 
assumed to be equal to unity and the velocity ve was calculated. 
Homeaneexact solution ve should be equal to unity assuming 
that the velocities along all the free surfaces are taken 
Squall tO carey: Evidently, We Wiebe GUI “OutY TO be: Larger 
than unity Beceee of the fact that the pressure along AB is 
larger than zero and thus the fluid accelerates along BC and 
decelerates along AL. The velocity at C was immediately 
€errectea to a first order of approximation by iterating on 


ve through the use of the following FORTRAN STATEMENT: 


V@ = V@ - (VC - 1.0) * V@ Crake: 


fete 
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iieemerocess Of iteration has been used in all subsequent 
imperati1ons. Apparently, the iteration on NP could have been 
eccelerated by beginning the analysis with a We value smaller 
mmen unity (say, 8 = 0.9). 

subsequently, the velocities along the free boundaries 
were calculated everywhere through the use of the $ values 
eeeeme NOdal points on the free surface and, as anticipated, 
it was found that they were not everywhere equal to unity. 
The reduction of these velocities to unity within a prescribed 
error constitutes the essence of the iteration process and 
tne use of the moving grid system. 

Let ve mepresemueune velocity at a point alone the 
teem suriaces. Then, the free surface at that point was 
moved inward or outward, depending on whether Le was greater 


wimeomaller than unity, by an amount given by the following 


statement: 
DELV = (VS ** 3 - 1.0) * FAC Cy) 


in which FAC is a Mw C up later Ine correction DELVewas appl red 
vo che ace POinvcme (Onpwemlime Nornal- or nearly normal vo 
the boundary. Various values of FAC have been used in the 
pretmanary caleulations, and it was found that a FAC value 
smaller than 0.05 will yield a convergent iteration. In the 
present analysis either one of the following two techniques 
have been used in correcting the free surface. In the first 


method a single FAC multiplier of 0.015 was applied to both 
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the upper and the lower boundary, five iterations were per- 
fermead, ana then the angle of deflection 8§ was increased or 
decreased by two degrees depending on whether the velocities 
Sia tne upper free surface was smaller or larger than unity. 
This process has been continued until everywhere the veloci- 
fms reached avvallue Ofel t+ 0.01. It is important to note 
that the angie of deflection, the jet thickness, and the 
rmadial coordinate rp, ond rq change qdurime each ei temacion pri- 
marily because of the corrections made on the free surfaces. 
mowever, the change in gp per iteration was about 0.25 degrees. 
Consequently, one would have needed at least 40 iterations to 
reach a correct jet position from an initially assumed posi- 
feoneweon an error of + 10°. Evidently, then, the + 2° 
increments applied to 8 have accelerated the convergence 
Process. 


In the second method two FAC values were used; namely, 


iH 


FACU = 0.001 and FACL We lee oven LLerac1OnSs were icCarriee 
out by applying FACU to the upper portion of the jet (along 
GH in Fig. 14) and FACL was applied to the lower boundary. 

At the end of the fifth iteration the role of the FAC values 
were Ae onaneed ands hivemmeCres.Feraui ons were carried OuL. 
Tats simply amounted to first iterating on the lower boundary 
manure n On ewute Upper boundary. As before, the deflection 
angle was incremented by two degrees until the correct jet 
deflection angle has been arrived at, as evidenced by the 


uniformity of the velocities along the upper and lower free 


surfaces. 
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Micwiwuimeescrintion Of the iteration process would 
not be complete without stressing few other computational 
Mme@ecalres which have enabled the prediction of the initial 
Gen tection angle and signaled at the arrival of the correct 
jet position. These will now be described. 

ine boundaries were iterated upon through the use 
memeaciier One Of the two methods cited above until a deflec- 
tion angle has been arrived at which was not more than a few 
@eerecs Off from the final jet position. At that time no 
limitations were imposed on Cae sump er Of iterations and the 
wereevas allowed to arrive at its final position. In passing, 
ipeeomould be noved that the programs were written in FORTRAN 
iyeana double precision arithmetic was used. Since the pro- 
grams on the Naval Postgraduate School computer (IBM 360-67) 
nedee Oo De run Overnight for a period of about one hour, the 
stops between two successive runs have provided the writer 
Weecimeeat Opportunity to examine the variation of the veloci- 
tiles, to observe the speed of convergence and finally to 
eeecide as to whether one has arrived at the stage of letting 
ie computer finely carry Out twenty to twenty-five itera= 
menos tO Pee abl i sti theo rina jel. pOSsit10N. 

Suwdugari iy. in dealing with the analysis of jet de= 
flection problems from arbitrarily~shaped curved surfaces, 
one will have at first very little insight as to the value 
Stmtncmactlcer1on ancie which the jet wall finally assume. 
However, after one or two calculations with one or two geo- 


Mevmuca ly combinations Of a given bucket and nozzle shape, 
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Piemva  enov Only be able to guess fairly accurately the 
Hepat jeu-Geflection angle but will also be able to derive 
a predictor equation. We will now describe the derivation 
of such an approximate equation which has been useful in 
preeaicking an initial deflection angle for all the remaining 
peses POLlOwing the study of the first case. 

Comsrcer the eduatvion of momentum applied to the 
fluid domain under consideration: 


2 2 2 2 _ 
a aac. 2m) pV,V,cos B= T K13) 


in which P, represents the uniform pressure along AB, V; the 
velocity along the free streamlines ole = an vane wpresene 
analysis), pop the density of fluid, and T the force required 
memaola tie bucket at rest. Evidently, if Mise B, and we 
were Known, then the force acting on the bucket could have 
been calculated. In fact, the present analysis does enable 
One wero dO SO not only through the use of the momentum equa-— 
MO onmotisor cheouen the dntearation of P dA. (aA 28-the 
projected elemental area of the bucket at the center of which 
the pressure is P). Since the velocities along the bucket 
were known for the first case, P dA was integrated and the 
resulting thrust was calculated. This procedure has not 
only provided a means of checking the thrust calculated from 
the momentum equation, but also an equation to predicu (EMe 
deflection angle for other geometrical combinations of the 


HOgeVemande rine. buUeCKeL., “nus, writing 
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V V 


faraAk = 2 roe e + (G2)® + a cos 8B C7) 
J J 
one has 
iL O 
fg f PdA - wll + Nae j 
cos B = — (75) 


2m(V/V 5) 


The equations so far derived are exact and Eq. (75) could be 
used as a mean to predict the angle of deflection for other 
ee and Sy Jin values if one assumes that f PdA remains 
essentially constant. Thus using the f PdA value obtained 
meen vie first case and a ve value of approximately 0.8, 

one arrives at the initial value of 8. It should be stressed 
Mniavyv the foregoing procedure is designed simply to help the 
fee se lCO ofr ye Taster at the correct solution. It neither 
alters the final result nor is absolutely necessary for the 


solution of the problems discussed here. 


D. HEMISPHERICAL TARGET-TYPE THRUST REVERSERS 

The use of thrust reversers on modern jet craft for the 
purposes of decreasing the landing load, maintaining 100% 
engine speed while landing, decreasing the thrust while 
maneuvering, reducing the length of the landing strip, tire 
wear and tear, etc. has been one of the most important 
developments in aircraft industry in the last two decades. 
To perform effectively, a thrust reverser should yield a 


relatively small jet-deflection angle, should have easy 
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Beowage configuration, should have relatively small boattail 
or stowage drag, and finally should not adversely affect the 
eeeine performance. A hemispherical thrust reverser appears 
to fulfill these requirements. Furthermore, there exist 
some experimental data on this type of reverser (Steffen et 
Pewee 34 | and Povolny et al. [35]) for comparison with the 
feeoults Of Che present analysis. 

»5o far as the finite element method is concerned, there 
Mmmeceresericu2on on the nozzle and the bucket shape and they 
eeuld be of any axisymmetric, reasonably-chosen configurations. 
Mievner words, the finite-element method enables one to 
feos Ze, Within the assumptions of the inviscid flow theory, 
Mies relavive position of a given bucket and nozzle, the effect 
of some geometrical modifications on the bucket and/or nozzle 
shape (such as placing a flat plate within the bucket) on the 
Permiermance of the reverser. Furthermore, it enables one to 
assess the degree of significance of the real-fluid effects 
such as viscous and turbulent energy dissipation, entrainment, 
Coanda effect, nozzle-pressure ratio, etc. on the actual 
performance of the mieere 

The geometrical configurations, namely, the ratios Ria alee 
and S/T. were chosen partly on the basis of the past experi- 
ence with such reversers. For example, it has been experi- 
mentally demonstrated by Steffen et al. [34] that the larger 
bucket, the larger is the reversed thrust and that the closer 
the nozzle to the bucket, the larger is the back-pressuring 


effect. In other words, when a nozzle is placed very close 
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wo Utne bucket, the pressure rises in the nozzle because of 
iemolockage Of the flow by the bucket and this, in turn, 
fonds tO aitect the engine performance. Therefore, there is 
Pemeprimum nozzle—-bucket configuration which yields optimum 
feverser performance. 

iis On the basis Of such considerations that the nozzle- 
Ducket feometries listed in Table I have been chosen for the 


numerical analysis and experiments. 


TABLE I 


Geometrical Charact. Angle of Deflection B (deg.) 


of the Hemispherical 


Buaiee ead Exoerimental 
ee aS Hance velement analysis (water) 
slope of mean slope of slope of 
the inner deflection the outer the outer 
Case R/T, s/P surface angle surface surface 
No. 
all ale: On4 fae t Z0e3 Ore) alk 
2 alpee O26 Shed Bou 0 Boo 50 
3 lege: ers Soe S055 Hi. 6 6 


The coordinates of the free surfaces for the three cases 
tive migmehaate [ware presenved in Appendix H. The computer 
program used in the study of the case number 3 is given in 
Appendix I. Suffice it to say that similar programs were used 
fOrwtmemanal veils Of the Olher two cases. The minor differences 
between the three computer programs essentially resulted from 


thesuse of different moving and fixed grids to suit the 


84 





Poeerwcular peometry Under study. An effort has been made to 
er euicemunper sol Clements and the bandwidth in the resulting 
matrix as low as possible and yet maintain a reasonable 
degree of accuracy. At various stages in the evolution of 
ie "COMpmmet sor OPramc. larger or smalller number of elements 
have been used, and it was found that the number of elements 
used herein (288) is quite sufficient. Needless to say, the 
coordinates of theefeee sumfaces may be calculated to a 
greater degree of precision by increasing the number of the 
EwenenivucweareLeularly Imetne moving ¢rid. The additional 
eorrections to be applied to the nodal points on the free 
Sur iacCameeumenMecncheNol “Unewitmial 1teratron were less than 
107°, It is estimated that the calculated deflection angle 
Hill notevercwemore than £1,0°. Sim view of the fact that the 
real-fluid efitects such as entrainment, Coanda, etc., can 
decrease the deflection angle by as much as ten degrees, it 
was not deemed necessary to pursue the calculation of the 
G@erleevion angle to a greater degree of precision. 

A series of experiments were conducted with water through 
the use of aretient. sharp-edged nozzles of appropriate 
Gemensions (see Table 1) with a hemispherical cup of 0.72 
inen ~radius primarily for the purpose of determining the 
Slope of the outer surface of the deflected jet. The slope 
was determined from photographs similar to that shown in Fig. 
ES end tanwlaved, in column {of Table. J. 

It is evident from Table I that the experimentally detver= 


Tmcdevalies or Bare Somewhat Smaller than those obtained 
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Figure !5 water jet deflecteg by a nemispherical bucket 


(Roy ete sl, Oy S/ 6.) = 2058) 
OO O 


es 
Cc 





memerically (see columns 6 and 7, Table I). The reasons for 
these differences are several, the most notable ones being 
Meewcnergy dissipation due to viscous action and the Coanda 
meececo. ine Jatvver will be discussed in detail parener later 
mmeeoniect10On with the experiments conducted with air. 
Suffice it to say that the Coanda effect with water is 
relatively small and the difference between the computed and 
the measured angles of deflection is within the range of 
experimental errors and numerical AOoCo eG ome 

The theoretical and experimental results may also be 
Compared in terms of a measure of the reversed-thrust because 
@oeaus ObVious practical significance. 

iaeone aircrait industry the performance of the reverser 


is expressed in terms of a reverse-thrust-ratio Np defined by 


 aetual reversed Thrust 


os 


Np forward thrust of the nozzle alone 


Cin Other words by 


The Np values found in the present analysis, and those 
obtained experimentally by Steffen et al. [34] are tabulated 


below. 
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TABLE II 








Case n n 

No. HM x “ 
Present Analysis Experimental [34] 

1 gy fase TOs 

2 O67 25 Oar de, 

5 Geese, 0.62 


Evidently the Np values obtained experimentally by Steffen et 
al. [34] are somewhat larger than those obtained in the present 
analysis. As noted earlier, the observed differences are 
primarily due to the Coanda effect and the nozzle-pressure 
ratio. The Coanda effect decreases the deflection angle and 
miemepy increases the reversed thrust. In the experiments 
conducted by Steffen et al. [34] the outer surface of the 
mozzZle was streamlined in the form of a boattail to decrease 
entrainment and te WO nde Credgewp wOumh OL ncreaac Np: In 
fact, as noted by Steffen et al. [34], the pressure reduc- 
PrLOns ON ane boattail were large enough to account for as 

much as 10 to 15 degrees reduction in the deflection angle. 
Thus the relatively large differences between the Np values 
computed in the present analysis and those obtained by Steffen 


et al. [34] are primarily attributable to the Coanda effect. 





MimouGcweco aiilustrave the influence of the Coanda effect, 
See particular experiment was carried out with air using a 
iee> inch radius nozgzle and a 2 inch radius hemispherical 
Picket with a 1 inch spacing (i.e. with Ae =m. and 
2/ = 0.8). The photograph of the deflected jet (visualized 
with aluminum powder) is shown in Fig. l6a. Figure 16b shows 
the deflection of a water jet from a similar nozzle-bucket 
geometry. It is evident from the comparison of these two 
photographs that whereas the slope of the outer free surface 
with air is 33 degrees, with water it is 46 degrees. The 
13 degree reduction in the deflection angle due to the Coanda 
effect corresponds to an approximately 20% reduction in Np 
Mims Drings the calculated Np values into closer agreement 
with those experimentally obtained. 

The nozzle-pressure ratio (total nozzle pressure/ambient 
pressure) or the actual velocity in the nozzle causes varia- 
tions in 6 primarily because the entrainment needs of the 
deflected jet and hence the Coanda effect increases with 
increasing jet velocity. 

SUL hace? 1c) CO say Vieveune analysis Ol thrush reversers 
is extremely complex for the cases where the jet impinges 
asymmetrically or where the gravitational effects are impor- 
tant. Furthermore, some of the practical problems associated 
Wit thrust reversal, suweh as the reattachment of the jet vo 
the nacelle of the engine, hot gas reingestion, interaction 
of the deflected jet with the ambient stream and the adjacent 


surfaces (in this case the outer surface of the nozzle) cannot 
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(b) Water jet 





Figure 16 Air and water jets defiected by a hemispherical 


bucket Layton leo, S/T. =80.c) 
0 O O 
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yet be analyzed by any one of the existing numerical methods, 
even for relatively simple two-dimensional geometries 
Composed of straight boundaries. Entrainment depends on the 
Menocity distribution and the turbulence level in the jet, 
igpe Ol the adjacent surfaces, the physical properties of 
the jet and of the fluid medium into which the jet discharges. 
Most of these parameters are obviously interdependent and at 
meilesent there does not seem to be any hope that their effect 
Smee ageflection of an axisymmetric jet can be evaluated. 
Nevertheless, approximate analyses such as the one 
presented herein help to isolate the more promising types of 
Maeust reversers, investigate their potentialities, evaluate 
mem ideal periormance characteristics, and to delineate 
Miemranete Of importance of the geometrical variables involved. 
Then those cases that are shown to be promising through 
Gomputer experiments can be tested in a laboratory, at a 
sufficiently reduced expense, to evaluate the effect of 
additional fluid-mechanical phenomena on their actual 


wert Oor Malice. 
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IV. SUMMARY OF RESULTS AND CONCLUSIONS 


Pectass Ol “allalytical and numerical solutions to the 
problem of inviscid-jet deflection from plane and axisym- 
metric concave surfaces have been presented. The methods 
Of solution for the two-dimensional cases included the 
modified-hodograph method developed herein as well as those 
devised by Levi-Civita and Riemann and Hilbert. 

The modified—-hodograph method is capable, within the 
Jimitations of the potential~flow theory, of analyzing the 
frow in curved obstacles and generating a family of practi- 
cally suitable bucket shapes. Furthermore, this method has 
Wem mrnmaerehny capability of being extended to the generation 
of even a larger family of concave-obstacle shapes through 
icine WO1 SGletvrapuced SOurces and Sinks in the nodograph-= 
plane. 

The method devised by Levi-Civita and successfully 
extended and employed herein requires the judicious selection 
Ot the variation of two parameters and the use of a suitable 
series (such as the Fourier series) to enable one to perform 
the integrais leading to the coordinates of the obstacle as 
well as of the free surfaces. A brief familiarity with this 
method allows one not only to generate a family of suitable 
ConGave Obstacles butnalso to round the corners of a buckere 
composed of straight segments. This procedure has been 
illustrated by rounding the corners of a U-shaped obstacle 


previously analyzed by Sarpkaya [1]. 
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As to the Riemann—Hilbert mehtod, it has been demonstrated 
biiec it 1s in essence a mathematically more elegant and 
sophisticated version of Levi-Civita's method in prescribing 
SereoiierUncer1ons. It is for this reason that no particular 
solutions have been presented and only a linear variation of 
the parametric function was discussed. 

These methods have provided indirect solutions for the 
Sliape Of the curved solid boundary in terms of a given jet- 
@emorcure and jet—deflection angle. The direct solution of 
the corresponding problems requires the use of either a 
Sequence Of 10eracions or a numerical scheme for solving the 
Maptace-Iield problem. The latter method has been reserved 
i@emumenanalysis of the axisymmetric cases because of their 
mememreaching and relative practical signiticance. 

The axisymmetric cases, namely, the deflection of a jet 
fmcudne froma Nozzle from hemispherical thrust reversers 
mavembeen Solved directly through the use of the finite 
element method. In doing so, several new procedures have 
been introduced. These include the use of a self-adjusting 
anda fixed grid and a convergent iteration scheme. The jet- 
deflection angle, the slope of the inner and outer surfaces 
of the deflected jet as well as the coordinates of the free 
streamlines were determined for a series of representative 
nozzle-hemisphere combinations. It has been demonstrated 
that the method and the iteration scheme are sufficiently 
Cenemwa hk tOsanalyae tne More promising Gyoes Of Thrust—-reverserc, 


evaluate their performance characteristics, and to delineate 
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the range of importance ofthe geometrical variables involved. 
Finally, a few experiments have been carried out with air 

ecu Waber FCO assess the degree of significance of the real- 
dived effects on jet deflection. These experiments as well 

as those conducted by others showed that in cases where the 
Coanda effect is negligible (water jet in air), the analytically- 
predicted and experimentally-measured jet-deflection angles 
are in good agreement and that the real-fluid effects are 
Meeiteible. This is evidence of the fact that the deflection 
of a free jet by a solid boundary is well suited to potential 
flow analysis particularly for the cases where the inertia 
forces dominate the establishment of the flow pattern. In 
cases where the Coanda effect is not negligible (air jet in 
air), the real-fluid effects are more pronounced and can cause 
the deflection angle to be as much as 20% lower than those 
predicted from the potential-flow theory. Nevertheless, the 
potential-flow analysis of the two-dimensional as well as 
axisymmetric jet-deflection problems enables one to determine 
the optimum values of the geometrical variables involved and 
thereby eliminates the use of costly trial and error techniques. 
With such information available, the real-fluid effects could 
be evaluated with a few carefully chosen experiments and the 
refinements of design could be based upon a secure knowledge 


Of tne Moowm prancaples ol mydrodynamics . 
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APPENDIX C 


EVALUATION OF THE COEFFICIENTS OF THE 
FOURIER SERIES REPRESENTING 6(a) 


Let o(co) be represented by the series 
oO 
Q(o) = ¥ a cos no aieceet MO RAuilana ame ara (A-1) 


Mecording to Fig. 6, @(q) is defined by 











O0(0) = = (ie cance O41 
re) 
Bh 
e(o) = F fo 2 | —— <PGs-— 
2 05705 T5045 i 2 
O(c) = te 3 CO. Oa 5 (A~2) 
ae ui 2 
A(c) ac Be CS TW Aele 
@(a) = o. ie eee 1-650 0 -< 1-6 
€ |oo704 Oo" 04 2 1 
@(o) = 0 TO, <0< 7 


Multiplying Eq. (A-1) by cos no and integrating over a period 


ene has. 


TT T 


f @(0) cos no do = a. if oa no do (A-3) 
—T —T 


Now making use of the symmetry, one obtains 


6(co) cos no do | (A-4) 


jab) 

ii 
= |N 
C2 


oH 





Tne above equation can be written as the sum of six integrals 


corresponding to each range defined by Eq. (A-2). Thus, one 





has, 
2 vd 2 aT O al 
a=-—/f (-n)cos no dg t+ f = Se ea raee| |) ees Glen (66) 
me 1 0 o, o 555, 05 07 
1 
2 Lids 
+f (- 5) COs no do + f (— 5) COs no do 
Oo JE 
O 2 
panath , 1c 1 
eal) = eae eee cos no do + f QO cos node 
=o 2. A Zoe T-O 
Ps 1 
(A-5) 
Mieber the 1ntbegrataon one obtains, 
7 1 n 
a = ———[ (cos now - cos no_) - (-1)° + (cos now - cos no,) | 
C 2 1 Pa ih 
n (5-0) 
(A-6) 
DOR = pai 3 ners 
or 
es oe (cos no. - cos no,) 
ig| (65-64) (A-7) 


iO eee. == lei. 7 
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mer n= 0, one has 
—_—[-| —_S-_CO _— 
a 5 (A-8) 


as expected from the average value of the function @(0). 


og 
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APPENDIX E 


EVALUATION OF THE INTEGRAL FOR w(t) 


_ Netter) 8 _ -(atbn) 
i. al (n-t) y¥-n(n-k) 


w(t) dn (E-1) 


Newnes that 


a by sya one 
il gee or ee ate 
n= n-t n- Tt ° 


Eq. (E-1) can be written as, 


k k 
Cae 9) {(atbt) ¢ Th + b / i | 
: lL (n-t) f-n(n-k) i. Sane 


(eo) 


The second integral in Eq. (E-2) is 


Tees 8 eo (E-3) 


The first integral in Eq. (E-2) depends onthe particular range 


Otieger  MOriG 6 leon One Nas. 
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K ee 
7 an : 1 7 Coal —oG + Vk-1 Vit (k—t) 
1 @-t) fnGen) fe 0et) = (t-1) 
(E-4) 
Herey < 0 and t > k one obtains, 
k em - -t+1 
cos © = (E-5) 


dn al 
ee ee 
1 (n-t) yn(k-n) / t(k-t) t-1 


Inserting Eqs. (E-3), (E-4) and (E-5) in Eg. (E-2) one finally 


GWirains 5 


_ atbt 2 el =o + \k-1 \/t (k-t) by 


w(t) -—| leat) sae ones 
= £ (t - 1) 
hom O < tsk, and 
io, eal -— 
1 ANE | TI k 


Tor t <—{05ond. tee kK 
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APPENDIX F 


EVALUATION OF THE DEFLECTION ANGLE B FROM EQ. (52) 


Letting t+o- in Eq. (52) one has 


k 
lim w(t) =i sy fatbnd g, (F~1) 
tre 1 Vn(k-n) 
and noting that 
B = telim iw(t) (F—2) 
-?00 


Eq. (F-1) can be rewritten after some rearrangement as, 





bk k 
(a +—-=) k 
B= nt—S— ps SL _ 2 n(n) (F-3) 
bs Vesicle, 1 


Bvaluating the integrals and applying the limits of integration, 


one nas , 





| (a +B cost <a + pia | (FH 4 ) 


lo 
ul 
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